Pion and kaon valence-quark parton distribution functions 
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A rainbow-ladder truncation of QCD's Dyson-Schwinger equations, constrained by existing appli- 
cations to hadron physics, is employed to compute the valence-quark parton distribution functions 
of the pion and kaon. Comparison is made to tv-N Drell-Yan data for the pion's w-quark distribution 
and to Drell-Yan data for the ratio UK{x)/uTr{x): the environmental influence of this quantity is a 
parameter-free prediction, which agrees well with existing data. Our analysis unifies the computa- 
tion of distribution functions with that of numerous other properties of pseudoscalar mesons. 
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Experimental information on the quark and gluon par- 
ton distribution functions (PDFs) in the pion have pri- 
marily been inferred from the Drell-Yan reaction [IH5] 
in pion-nuclcon and pion-nucleus collisions. Kaon PDF 
data exists in the form of the ratio uk{x)/u-^{x) [llll]. 
While the nucleon PDFs are now fairly well determined, 
the pion and kaon PDFs remain poorly known. Reference 
[5] reviews both the experimental and theoretical status 
of nucleon and pion PDFs. Since the pion is central to 
hadron physics, and its key characteristics are dictated 
by dynamical chiral symmetry breaking, pion structure 
is a critical testing ground for our understanding of non- 
perturbative QCD. Much more theoretical work has been 
devoted to the pion elastic charge form factor (e.g., [B]); 
TTTT scattering (e.g., i?]); and the pion electromagnetic 
transition form factor (e.g., [SI) than has been devoted 
to the pion PDFs. Herein we take a material step toward 
ameliorating that deficit. 

Numerical simulations of lattice-regularized QCD are 
restricted to the computation of low-order moments of 
the PDFs: the pointwise dependence is not directly acces- 
sible O H] . Model calculations of PDFs are challenging 
because the framework employed must necessarily pos- 
sess perturbative QCD features in coexistence with the 
covariant, nonperturbative description of a bound state. 
Chiral symmetry has guided studies of pion PDFs within 
the Nambu-Jona-Lasinio (N JL) model [TUl [H] at the ex- 
pense of: an unphysical point-particle structure for the 
pion Bethe-Salpeter amplitude; and ambiguities from a 
dependence upon regularization procedure owing to the 
lack of renormalizability. Constituent quark models [12] 
and instanton-liquid models [13] have also been used. In 
all these approaches, it is difficult to have pQCD ele- 
ments join smoothly with nonperturbative aspects while 
respecting the quantum field theoretical nature of the 
underlying dynamics. The large x behavior of the pion 
PDFs provides an illustration. The QCD parton model 
[l4j and pQCD |15] are clear: at a scale of order-Aqco 
the behavior is Ut^{x) oc (1 — x)" with a = 2 + j where 



7 > is a logarithmic correction. However the above 
models imply an a ranging from to 1, or at most 1.5 [5]. 

These issues may in principle be addressed if the 
PDFs can be obtained from truncations of QCD's Dyson- 
Schwinger equations. The DSEs are a hierarchy of cou- 
pled integral equations for the Schwinger functions (n- 
point functions) of a theory. Bound-states appear as 
poles in the appropriate n-point functions; e.g., the 
bound-state Bethe-Salpeter equation (BSE) of field the- 
ory appears after taking residues in the inhomogeneous 
DSE for the appropriate color singlet vertex. Numer- 
ous reviews; e.g., |16j . describe the insight into hadron 
physics achieved through the use of the rainbow-ladder 
(RL) truncation of the DSEs, which is the leading-order 
in a systematic, symmetry-preserving scheme |17j . 

The first DSE study of PDFs was conducted for the 
pion [18" in an analysis that employed phenomenolog- 
ical parametrizations of both the Bethe-Salpeter ampli- 
tude and dressed-quark propagators. The purpose of this 
present work is, for the first time: to employ numeri- 
cal DSE solutions in the computation of the pion and 
kaon PDFs, adapting the RL model employed in success- 
ful predictions of electromagnetic form factors [6l fT9H2T] ; 
and study the ratio UKix)/uTr{x) in order to elucidate 
aspects of the influence of an hadronic environment. 

In the Bjorken-limit, DIS selects the most singular be- 
havior of a correlator of quark fields of the target with 
light-like and causal distance separation ~ 0+. With 
incident photon momentum along the negative 3-axis, the 
kinematics selects z"*" ~ z_l ~ leaving z~ as the finite 
distance conjugate to quark momentum component xP^ , 
where x — Q'^ /2P ■ q is the Bjorken variable, = —Q^ is 
the spacelike virtuality of the photon, and P is the target 
momentum. To leading order in the operator product ex- 
pansion, the target structure functions are proportional 
to the charge- weighted sum of PDFs, qf{x), for parton 



2 




FIG. 1. (Color online) Diagrammatic representation of par- 
ton distributions. Top panel - the exact parton distribution 
corresponding to Eq. (2]); and bottom - Rainbow-ladder trun- 
cation of the amplitude G. S2 is the qq propagator. 



of flavor /. The PDF is given by the correlator 
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expressed here in manifestly Lorentz-invariant form. 
In the infinite momentum frame, qf{x) is the proba- 
bility that a single /-parton has momentum fraction 
x = k-n/P-n |23]. In the above, n'^, and (for later 
use) p^, are light-like vectors satisfying = = and 
n ■ p = 2. They form a convenient basis for the lon- 
gitudinal sector of 4-vectors. One has k ■ n = and 
k ■ p — kr . The dominant component of q is parallel to n, 
i.e., q~ dominates. Note that qf{x) = —qjr^—x), and that 
the valence quark amplitude is qj{x) = qf{x) — qf{x). 

It follows from Eq. (jlj) that ^^dxq"^{x) — 
(7r(P)|J+(0)|7r(P))/2P+ = P^(0) = 1. Approximate 
treatments should at least preserve vector current 
conservation to automatically obtain the correct 
normalization. 

In our DSE framework, dynamical information on the 
various nonperturbative elements, such as propagators 
and bound state amplitudes, is available in a Euclidean 
momentum representation. (In our Euclidean metric: 
{it^.lA = 2(5^^,; -fl = 7^; 75 = 74717273; a-b = Y.j^^a.b^; 
fi = 7-n; and timelike < 0.) The corresponding 

formulation of Eq. ([ij is 



Cf^k 
(2^ 



5{k-n-xP-n) tv^A[ifiG{k,P)] 



(2) 

where trcd denotes a color and Dirac trace, and G{k, P) 
represents the forward q-target scattering amplitude. In 
Euclidean metric the vectors n, p, P satisfy n? — Q = p^ ^ 



n-p = —2, P'^ = —m^, and P ■ n 

The top part of Fig. [l] illustrates Eq. ^ . In rainbow- 
ladder truncation, which treats only the valence qq struc- 
ture of the target, we have the decomposition illustrated 



in the bottom part of Fig. [T] The 4-point function S2 is 
the qq two-body propagator and F^ is the Bethe-Salpeter 
bound-state amplitude, both computed in the RL trun- 
cation. The combination S2 ^ ifiS^k ■ n — xP ■ n) = 
S{e) r"{i; x) S{e), where F"(£; x) is a dressed vector ver- 
tex. The RL truncation for the valence UTr{x) is thus 



u^ix) = - 



2 J {277)4 

X Suii) F"(^; x) s^{e) r^{e, p) Sd{i - p)] , (3) 

wherein F"(.^; x) is a generalization of the dressed-quark- 
photon vertex, describing a dressed-quark scattering 
from a zero momentum photon. It satisfies the usual 
inhomogeneous BSE (here with a RL kernel) except that 
the inhomogeneous term is ifiS{£-n—xP-n). The dressed- 
quark propagator is 5(^;C) = l/[ifAie^;0 + B{e;C)], 
where C is the renormalization mass scale. 

This selection of RL dynamics parallels precisely the 
symmetry-preserving treatment of the pion charge form 
factor at = 0, wherein the vector current is conserved 
by use of ladder dynamics at all three vertices and rain- 
bow dynamics for all three quark propagators [51 
Equation ^ ensures dx q"f{x) = 1 for f = u,d auto- 
matically since / dxF"'(^;x) gives the Ward-identity ver- 
tex and the result follows from canonical normalization 
of the BS amplitude. 

We adopt the representation = ^{ap^ + /3n'^) + k^ 
to transform to new variables a = —£ ■ n and (3 = ~£ ■ p, 
thus converting Eq. ^ to the form 



u^{x) = 
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. (4) 



where: Je = —i/2 is the Jacobian of the variable trans- 
formation; and T is the result of the trace in Eq. (|3|, 
using F"(£; x) « nf,dS-\£)/dif, S{i-n-xP- n), which 
is the correct result from the Ward Identity after / dx. 

Since qf{x) is obtained from the hadron tensor W'^'^ , 
which in turn can be formulated from the discontinu- 
ity T'^'ie) - r'"'(-e) of the forward Compton amplitude 
Tf^"^ , we observe that all enclosed singularities from the 
difference of contours cancel except for the cut that pro- 
duced the delta function constraint on a. 

We employ the RL-DSE model developed in Refs. [24]- 
[26] . in which the BSE kernel takes the form K = 
—Att acs{k'^) D^i^°°{k)^j^ <Si ^Jv, where k is the gluon 
momentum. Here acff(fc^) is a model running cou- 
pling chosen such that it reproduces QCD's one-loop 
renormalization-group behavior for k^ >2 GeV^ . A more 
general method for treating K has recently become avail- 
able [27]. The DSE that produces the dressed quark 
propagator is also determined by (fc^) |24n26j : and the 
combination of the DSE and BSE produces dressed color- 
singlet vector and axial-vector vertices satisfying their 
respective Ward-Takahashi identities. This ensures that 
the chiral-limit ground-state pseudoscalar bound-states 
are the massless Goldstone bosons from dynamical chiral 
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TABLE I. Illustrative selection of DSE results P [T9l [20l [26] 
obtained with the RL kernel employed herein compared with 
experimental values 29 . (Dimensioned quantities are listed 
in GeV or fm^, as appropriate.) 
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expt. 
calc. 


0.138 0.092 
0.138 0.092 


0.496 0.113 
0.497 0.110 


0.44 0.34 
0.45 0.38 


0.5 0.42 
0.5 0.41 



symmetry breaking [231 US] ! and it ensures electromag- 
netic current conservation |28j . This kernel is found to be 
successful for, amongst other things, light-quark meson 
properties f2B] including electromagnetic elastic [BJ [TS] 
and transition HI] form factors. The model param- 
eters are the two current quark masses and one infrared 
strength for aes{k^)- Selected results related to the pion 
and kaon are displayed in Table |T| 

In the evaluation of Eq. Q we employ the full pseu- 
doscalar meson Bethe-Salpeter amplitude 

TA£, P) = To [iEAq; p) + PFAq; p) 

+ iGAq\P) + (T^,q^P,HAq\P)] , (5) 

where q = £ — P/2 is the relative qq momentum appropri- 
ate to Eq. ([3|. For a charge-conjugation eigenstate (e.g., 
the pion), the invariant amplitudes E, F and H are even 
in q ■ P, while G is odd. The kaon invariant amplitudes 
contain both even and odd components. We expand the 
q ■ P dependence in Chebschev polynomials 126' , keeping 
terms of order n = — 3. The domain of i'^ over which 
the quark propagators are needed in this application is 
larger than what is available from previous solutions of 
the quark DSE. We therefore adopt a constituent mass 
pole approximation for the denominator of the specta- 
tor quark propagator [TB]. A constituent spectator mass 
of 0.4 GeV permits a minimal adjustment to establish 
the normalization (x"). We compared the approximation 
r"(^; x) w nf,dS-\e)/d£f, S{£-n-xP- n) with the bare 
vertex truncation and found that both give essentially 
the same distributions after re-enforcing the normaliza- 
tion {x°) = 1. 

In Fig. |2] we display our DSE result [3D] for the va- 
lence w-quark distribution evolved to = (5.2 GeV)^ 
in comparison with ttN Drell-Yan data [3] at a scale 
^ (4.05 GeV)^ obtained via a LO analysis. Our dis- 
tribution at the model scale Qo is evolved using leading- 
order DGLAP. The model scale is fixed to Qo = 0.57 GeV 
by matching the x" moments for n = 1, 2, 3 to the exper- 
imental analysis given at (2GeV)^ [53]. Our momentum 
sum rule result 2 (x) = 0.74 (pion), 0.76 (kaon) at Qq 
shows clearly the implicit inclusion of gluons as a dy- 
namical entity in a true covariant bound-state approach. 

In Fig.|2] we also show the result from the first 
DSE study jiSj, which employed phenomenological 
parametrizations of the nonpcrturbative elements. Our 
present calculation lies marginally closer to the Drell-Yan 
data in Ref. [3^ at high- a;. However, this is not significant 
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FIG. 2. (Color online) Pion valence quark distribution func- 
tion evolved to (5.2 GeV)^. Solid curve - full DSE calculation 
[30] : dot-dashed curve - semi-phenomenological DSE-based 
calculation in Ref. [TH] ; filled circles - experimental data from 
Ref. [3], at scale (4.05 GeV)^; dashed curve - NLO re-analysis 
of the experimental data ^31 j and dot- dot- dashed curve - 
NLO reanalysis of experimental data with inclusion of soft- 
gluon resummation [32] . 

because both DSE results agree with pQCD; viz., u{x) ^ 
(1 — x)" with a > 2 and growing with increasing scale, 
which is not true of the reported Drell-Yan data. 

Motivated by this, a NLO reanalysis of the data was 
performed [3IJ; and we also show that result at = 
(5.2 GeV)^ in Fig. [2] It does clearly reduce the extracted 
PDF at high-j: but not enough to resolve the data's ap- 
parent discrepency with pQCD behavior, which is dis- 
cussed at length in Ref. 0. The DSE exponents are 2.4 
at model scale Qq = 0.54 GeV in Ref. [T^, and 2.1 at scale 
Qo = 0.57 GeV for the present study. DSE analyses do 
not allow much room for a larger PDF at high- a; . A res- 
olution of the conflict between data and well-constrained 
theory has recently been proposed: a reanalysis of the 
original data at NLO with a resummation of soft giuon 
processes [32] produces a PDF whose behavior for x > 0.4 
is essentially identical to that of the earlier DSE calcula- 
tion [18], as is apparent in Fig.[2j 

In Fig. [3] we display the first nine moments of our re- 
sult for Ut^{x) at scale Q^ = (5.2 GeV)^ in comparison 
with the earlier DSE result from Ref. [IH] and the NLO 
reanalysis '3T| of the original E615 data, all plotted as a 
%-deviation from the moments of the most recent anal- 
ysis of Ref. [32 . Considering that the high moments are 
small, e.g., (x^) ~ 0.003, the two DSE results are both 
equally well in accord with the recent analysis. 

The ratio uk/ut^ measures the effect of the local 
hadronic environment. In the kaon, the w-quark is part- 
nered with a heavier partner than in the pion and this 
should cause u{x) to peak at lower- a; in the kaon. Our 
DSE calculation [30] is shown in Fig. [4] along with avail- 
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FIG. 3. (Color online) Moments of the pion's valence u{x) 
at scale (5.2 GeV)^, shown as a % deviation from the re- 
cent (ASV) re-analysis [32] (NLO, with soft gluon resumma- 
tion) of the 1989 E615 TriV Drell-Yan data 0]. Filled cir- 
cles - present full DSE calculation [30]; filled squares - semi- 
phenomenological DSE-based calculation [T5]; and filled dia- 
monds ~ re-analysis (NLO, without soft gluon resummation) 
of the same Drell-Yan data [31] . 



FIG. 4. (Color online) DSE prediction for the ratio of u- 
quark distributions in the kaon and pion [5] [30]. The full 
Bethe-Salpeter amplitude produces the solid curve; the re- 
duced BSE vertex produces the dashed curve. The reduced 
amplitude retains only the invariants and amplitudes involv- 
ing pseudoscalar and axial vector Dirac matrices, and ignores 
dependence on the variable q ■ P. These are part of the reduc- 
tions that occur in a pointlike treatment of the pseudoscalar 
mesons. The experimental data is from [311]. 



able Drell-Yan data [U H] , which does not separate sea 
and valence quarks. Our parameter-free result agrees well 
with the data. It is notable that the value of this ra- 
tio at a; = approaches one under evolution owing to 
the increasingly large population of sea-quarks produced 
thereby [31] ■ On the other hand, the value at x = 1 
is a fixed-point under evolution; i.e., it is the same at 
every value of the resolving scale Q^, and is therefore a 
persistent probe of nonperturbative dynamics j5j. 

In Fig.|4] we also display a calculation which employs 
a reduced BSE vertex: only the leading two invariant 
amplitudes E{q\ P) and F{q\ P) are retained, and each 
is truncated to the lowest Chebychev moment in q • P, 
i.e., E{q;P) E{q^). The field theory variable q ■ P is 
a constant in quantum mechanics. (These reductions in 
the BSE vertices occur within a NJL model description; 
but that model also ignores the q"^ dependence of the 
vertices.) These simplifications do not change the quali- 
tative behavior of the ratio, but the detailed quantitative 
agreement is impaired. 

The ratio in Fig.|4] was also computed using a 
nonpointlike-pion-regularized version of an NJL model 
[lOj . In this case one finds that while the behaviors of 
u^ix) and u^{x) are separately poor, the ratio is never- 
theless in fair agreement with extant data: the result is 
similar to the dashed-curve in the figure. 

An estimate of the leading large-x behavior 
uk{x) ^ Ak — x)" can be made in the limit 
where the quark propagators are characterized by 
constituent masses M„,Ms and the vertex is taken to 
be ifi5{l ■ n — xP ■ n), preserving the Ward Identities. 



We also truncate to j^EKiq'^) = J^Nx/iq'^ + ^k) 
where q — £ — P/2. The quark mass dependence of 
and At^ will provide an estimate of Ui<-(l)/u7r(l). For 
X > 1/2 the pole in the spectator quark propagator is 
the only one in the upper half plane and the i~ integral 
may readily be evaluated to yield 



UKix) 



dfi 



xM"' 



M, 



[M + M2]2 



- Eliq')- 



(6) 

Here: = m\ — {Mg — M^)^; we have changed the 
integration variable from to ji = ^ where the 
latter is the value at the pole; q^ evaluated at 
the l~ pole is = + {^J. — M'^)/2; and Hm{x) = 

a/ [1 — x) — xm?j^, with a = xM^. This divergence of 
the lower limit for large x guarantees that the result is 
completely determined by the ultraviolet behavior of the 
propagators and bound state amplitudes. 
The integral can be expressed as 



ux{x) — N 



dfi 



^1 



d + jiy 



(7) 

where bound-state amplitudes determined by one gluon 
exchange correspond to n = 2. The quantities a, 6, c, d 
depend on the mass-dimensioned scales in the system 
and are nonsingular in a;: a scales with the square of 
the spectator quark mass and other details are imma- 
terial. A change of variable to pL — [1 — x)fl/a shows 
that uk{x) oc [(1 — a;)/a]" when a/{l — x) is greater than 
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any physical mass-scale in the system. Running of the 
struck quark mass over a wide domain can be accom- 
modated. We thus have uk{x) oc Nx (1 — x)"^ /M^ and 
M^(a;) (X TVjr (1 — x)^/M*. Note that it is the bound-state 
amplitudes that completely determine the exponent a [5] : 
if the argument of Ek/tt did not diverge at large-x, the 
combined scaling effect of the propagators would vanish, 
giving a = 0. 

The above analysis applied to the ratio leads to the ap- 
proximate formula Ui<-(l)/u7r(l) ^ [Mu/Ms)'^ ~ 0.2, 
where the ratio of Bethe-Salpeter amplitude normaliza- 
tion constants is estimated from the experimental f-rrl fx- 
This estimate is in fair accord with our full calculation 
in Fig.|4] The NJL model with a sharp cutoff yields 
[Mu/MsY' [IQ]. However, in general this lacks a physical 
contribution from bound state amplitudes and NJL re- 
sults depend sensitively upon the regularization scheme. 

With this study we have unified the computation of 
distribution functions that arise in analyses of deep in- 
elastic scattering with that of numerous other properties 



of pseudoscalar mesons, including meson-meson scatter- 
ing and the successful prediction of electromagnetic elas- 
tic and transition form factors. Our results confirm the 
large- a; behavior of distribution functions predicted by 
the QCD parton model; provide a good account of the 
■n-N Drell-Yan data for Ut;{x)] and our parameter-free 
prediction for the ratio uii{x)/uT;{x) agrees with extant 
data, showing a strong environment-dependence of the 
M-quark distribution. 
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